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ABSTRACT 

This report attempts to correct some misunderstandings that have appeared in the literature 
concerning the order of accuracy of the QUICK scheme for steady- state convective 
modelling. Other related convection-diffusion schemes are also considered. The original 
one-dimensional QUICK scheme written in terms of nodal-point values of the convected 
variable (with a 1/8-factor multiplying the “curvature” term) is indeed a third-order repre- 
sentation of the finite-volume formulation of the convection operator average across the 
control volume, written naturally in flux-difference form. An alternative single-point upwind 
difference scheme (SPUDS) using node values (with a 1/6-factor) is a third-order 
representation of the finite-difference single-point formulation; this can be written in a 
psewdo-flux-difference form. These are both third-order convection schemes; however, the 
QUICK finite-volume convection operator is 33% more accurate than the single-point 
implementation of SPUDS. Another finite-volume scheme, writing convective fluxes in 
terms of cell-average values, requires a 1/6-factor for third-order accuracy. For complete- 
ness, one can also write a single-point formulation of the convective derivative in terms of 
cell averages, and then express this in pseudo-flux-difference form; for third-order accuracy, 
this requires a curvature factor of 5/24. Diffusion operators are also considered in both 
single-point and finite-volume formulations. Finite-volume formulations are found to be 
significantly more accurate. For example, classical second-order central differencing for the 
second derivative is exactly twice as accurate in a finite-volume formulation as it is in single- 
point. 


*Work funded under Space Act Agreement NCC 3-233. 
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INTRODUCTION 


Ever since the QUICK scheme 1 was introduced in 1979, there appears to have been some 
controversy regarding the formal order of accuracy of the convection terms. There have 
been attempts to clarify the situation 23 , but a recent journal article 4 indicates that 
considerable confusion remains in the CFD literature. This report is a further attempt to 
correct some of the misunderstandings that have arisen. 

First, for definiteness, consider a one-dimensional model problem on a uniform grid 
of mesh-size h, numerically simulating a (nondimensional) convection-diffusion equation 


d<j> 

dx 


1 d 2 <t> 
P € dx 2 


+ S(x) 


( 1 ) 


where Pe (= const > 0) is the macroscopic Peclet number. With appropriate boundary 
conditions, this represents steady-state one-dimensional convection and diffusion of a scalar, 
with a known source term. A finite-difference numerical approximation of this problem 
simulates 


(!£) = 

±| 


| \dxli 

Pe' 

idx 2 h * 


( 2 ) 


at each nodal grid-point, i. This will be called the single-point (SP) formulation. By 
contrast, a finite-volume formulation is obtained by integrating Equation (1) across a control- 
volume cell and dividing by h. This gives, for control-volume i, 



(3) 


where left and right face- values of the variable and its gradient are indicated, and S is the 
control-volume average of the source term. This is the operator-average (OA) formulation. 
Conservation of convective-plus-diffusive flux is guaranteed if the modelled face values and 
gradients satisfy 


and 


<M 0 

= 4> r d- 1) 

(4) 

*/<*) 

= VO'-D 

(5) 


Note that Equation (3) is an exact equation and not an approximation of Equation (2). They 
are related by (the one-dimensional form of) Gauss’ Divergence Theorem. 

In order to create numerical algorithms for (approximately) solving Equations (2) or 
(3), one needs to estimate either the derivatives in Equation (2) or, alternatively, the face 
values and gradients appearing in Equation (3). Additionally, in makijig these estimates, one 
has the choice of using either node- values, or cell-averages, <j> r For the operator- 
average finite-volume equation, the formulation will automatically be in flux-difference 
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form; numerical models of single-point formulations of derivatives can usually be written as 
the difference of terms satisfying conservation, thereby generating a pjewdo-flux-difference 
construction. 

The one-dimensional QUICK scheme is based on estimating face values and gradients 
using quadratic upstream interpolation through no</e-values of 4> located at the centre of 
control-volume cells. For example, at the right face, a parabola is interpolated through 4> i+ j, 
4>i, and (for P € > 0) giving the original 1 “1/8-factor” face value 

M r ) QU,CK = ■!«>,,, ♦ *,) - £(*,., - 2*. * (6) 


and, for the gradient, 


(</>/) QUICK 


iKx ~ *i) 

h 


(7) 


with left- face quantities obtained by lowering all indices by 1. Substitution of the QUICK 
formulae into Equation (3) gives, using classical Taylor-series analysis, 




0(h 4 ) 


1 

[*; - 

Pe 

h 


+ +«> 

24 ' 


0(h 4 ) 


+ S. 


( 8 ) 


In other words, in this finite-volume formulation, the QUICK convection operator average 
terms are third-order accurate, whereas the QUICK diffusion operator average terms are only 
second-order accurate. Thus, the overall QUICK scheme for the convection-diffusion 
equation gives OQi 2 ) convergence as the grid is refined. Controversy in the literature 
concerning the third-order convection term apparently stems from confusing finite-volume 
and finite-difference formulations. In particular, Equation (3) has sometimes been construed 4 
as an approximation of Equation (2). It is not. It is a perfectly valid (control-volume) 
equation in its own right. 


Operator-average or single-point formulations 

In order to clarify the distinction between finite-volume (operator average) and finite- 
difference (single point) formulations, consider the Taylor-series formulae 



II 

/ s 

4>, + X + i */' X 2 + 1 4>r X 3 + ... 


(9) 


4>'(x) = 

4>; ■+ 4>r X + i */" X 2 + I X 3 + ... 


(10) 

* * i *r 8 * 2 ♦ i *!"■” ♦ ... 

Now compute finite-volume formulae by subtracting the Taylor-expansions 
x = h!2 and x = -hi 2, giving 

(11) 
written for 


4>? ~ 4>f 

h 

.(„*!) + k 2 + /l 4 ^(n.5) + ft 6 ^0.*7> + 

‘ 2 2 3! 2 4 5! * 2 6 7! 

... 

(12) 
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This formula is also valid for negative n (representing integration); in particular, for n = -1, 
the control-volume cell-average of the transported scalar itself is given by 


+t 




*. + *! V 
V ' 24 


1920 


*? v) - 


(13) 


The expression on the left of Equation (12) represents the control-volume operator average 
(OA) of the (n+l)th derivative, whereas the first term on the right is the single-point (SP) 
form. Note that the difference between the two always involves an 0(h 2 ) quantity. This is 
an important point that is the key to clarifying the confusion that has arisen in the literature. 
If a finite-volume (OA) discrete operator is viewed as an finite-difference (SP) term, there 
is an 0(h 2 ) discrepancy between the two. This does not affect the leading truncation error 
of first-order schemes. Second-order schemes show a change in the numerical value of the 
/^-coefficient. But a third (or higher) order OA scheme is only second-order accurate when 
viewed as an SP scheme, and vice versa. This is apparently why the QUICK scheme has 
been so controversial. 


THIRD-ORDER-ACCURATE STEADY TRANSPORT 


In a recent paper 4 , Johnson and MacKinnon attempted to clarify the distinction between 
finite-difference and finite- volume formulations. Unfortunately, their conclusions are exactly 
the reverse of the true situation. They claim, in particular, that the QUICK(l/8) scheme is 
only a second-order accurate finite-volume convection scheme. For example (using here 
“left-right” rather than “east-west”), Johnson-and-MacKinnon’s Equation (4) for the 
QUICK(l/8) convection scheme is 

<i>r - 4>1 = 3 </> i+1 + 3 4>f - 74>j. 1 + 4>,-2 _ _ +."•) + HOT (14) 

h 8 h 16 ' 


This is correct as written; but Johnson and MacKinnon seem to imply, quoting Bradley et 
al. 5 , that this represents an 0(/i 2 )-accurate operator. First of all, from Equation (12), with 
n = 3, the leading truncation error is 



h'ltr’" ~ V" 

16 V h 


hi <t>™ + IL <t>p + 

16 384 


(15) 


Equation (14) is equivalent to the left side of Equation (8), showing the QUICK(l/8) 
convection term indeed to be a third - order accurate finit e-volume OA formulation of 
(4> r - 4> t )!h. Secondly, it must be stressed that the discrete operator in Equation (14) is not 
intended to be an SP numerical model of 4 >/ . If it is considered to be, as in Reference 4’s 
Equation (5), it will appear to be 0(1?) accurate; this is easily seen from Equation (12), with 
n = 0. 

Johnson and MacKinnon claim to demonstrate the “second ’’-order accuracy of the 
QUICK(l/8) convection terms by giving a numerical example of a simple convection-diffu- 
sion problem with a known exact solution, using a fourth-order accurate diffusion operator 
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„ = -_K2 + - 300,. * 16^,-t - <£,-2 + Qfh 4 ) (16) 

9i I2h 2 

the strategy being that the grid convergence will be dominated by the lower-order convection 
term. As is well known 6 , this is indeed a fourth-order finite-difference SP approximation of 
the second -derivative at point i, consistent with a quartic polynomial interpolated through 
node-points: <t> p 4> M , and </> i+2 . However, this is not what is being modelled in 

a flnit e-volume formulation. Rather, to be consistent, according to Equation (3), one should 
model the operator average across the control volume. The appropriate fourth-order finite- 
volume formula is 

4>r " = + 28 <t> i+l ~ 52 ^,- 28 ^,-l ~ ^i-2 + n(h 4^ (17) 

h 24 h 2 


More specifically, the fourth-order right-face gradient can be represented by 


tf>/(4th) 



^,♦2 " 3^*i + 30 f - 4> t _ x 

24 h 


(IS) 


obtained by interpolating a cubic polynomial through node-points: <£,, and 4> i+2 . 

The 0(h 2 ) convergence reported by Johnson and MacKinnon using QUICK(l/8) for 
convection and Equation (16) for diffusion occurs because their diffusion operator is only 
OQi 2 ) accurate in a finite-volume formulation. This can be seen immediately from Equation 
(12), written for n = 1: 







* hi 

24 ‘ 


— + 
1920 


(19) 


NUMERICAL EXAMPLE 

The numerical example used by Johnson and MacKinnon is (with a slight change in notation) 

d4> _ 1 d 2 4> _ 0 (20) 

dx P6 dx 2 

with boundary conditions on the nodal values 

m = o, <ki) = i (21) 

The exact solution is 

4>{x) 

Step sizes of 1/4, 1/8, 1/16, 1/32, and 1/64 are used, and Pd = 4. In the current 
formulation, pseudonode values are required beyond each end of the physical domain. For 
the purposes of this numerical test, these are taken here to be exact values given by 


, P«£jt 


- 1 


,P<S 


- 1 


( 22 ) 
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and 


+-i 


e- Ah - 1 
e 4 - 1 


(23) 


4 >»* i 


e 4(1 * w - 1 
e 4 - 1 


(24) 


Exacf derivatives and fluxes 

In order to investigate the effect of individual modelled terms, it is instructive to 
compute exact derivatives and fluxes from the known analytical solution. For example, 
errors introduced solely by modelled convection terms can be studied in isolation by using 
exact diffusion terms, and vice versa. From Equation (22), the first and second derivatives 
are 

d(j> _ P6 e Pf * (25) 

dx e ?i - 1 

and 

d 2 <j> _ Pe 2 (26) 

dx 2 e ri - I 


By integrating Equation (20) from (x - hi 2) to (x + hi 2) and dividing by h, the {exact) 
control- volume formulation is 



(<t> r ~ <t>t\ _ 

. 1 | 

->• 

1 

II 

o 

(27) 


l h I 

Pe' 

l h 1 


where 

4> t (x) 

= 

4> r (x - h) 

(28) 

and 

4> t \x) 

— 

<t> r ’(x - h) 

(29) 


By defining the convective-plus-diffusive flux at any point as 


F{x) 


<Kx) - <y(x)n>6 

h 


(30) 


Equation (27) can be written in flux-difference form across any control-volume cell of width 
h as 


F r (x) - F t (x) = 0 


( 31 ) 


where the exact convective-plus-diffusive fluxes are 
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and 


F r (x) 


F t (x) 


1 r {e Pix e PJ2 - 1) 

h e Pi - 1 


1 (e P6x e~ PJ2 - 1) 
h e Pi - 1 



e P6x e' p *' 2 " 
e pe - 1 . 


(32) 


(33) 


introducing the grid Peclet number, P A = h Pe. Note that conservation is guaranteed, since 
F t {x) = F r (x-h). 


Evaluation of discrete operators 

With a uniform grid of step-size h and a control volume centered at x„ a numerical 
model of the convective flux can be tested by using a hybrid formulation of Equation (30) 

*r“ - y] (34) 


and, assuming the numerical model to be conservative. 


I-.HYB 

F t 


(*,) 


77. HYB 
Ff 


(x r h) 


(35) 


If the modelled face values, tf 0 ** and <f>^\ are written in terms of nodal values of 4> 

<$>, = *(*,) ( 36 ) 

a solution of the flux-difference equation 

= 0 < 37 > 

then gives the computed <£, values corresponding to the particular convection model, treating 
the diffusive fluxes exactly . The node-point error is then, using Equation (22), 

NPE,. = 4> t (computed) - ( 3 ^) 


A grid-refinement study (with K 1 = 4, 8, 16, 32, and 64) then shows the true convergence 
rate of the convective model in isolation. The rate, R, can be obtained from 


| NPE.(A) | = 

|NPE 2i .(/i/2)| 

and should approach an asymptotic value as h becomes smaller and smaller. 


(39) 


I. Standard QUICK scheme 

Equation (6) for the right face value is rewritten here for convenience 


(* r ) QUICK 


HKx + *i) - t(*m " 2<t> ‘ + *1-1 > 


(40) 


with the left face value given by 
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0 M QUICK = T (*, + *m) - {(*, - 2 *.-i + *«) (41) 

Using exact diffusive fluxes as described above, the results in the QUICK(l/8) column of 
Table I show that the node-point error at x = 0.75 converges at a rate of 0(h 3 ). 


II. The SPUDS formulation 

An alternate formulation, recommended by some researchers 4 , is based on a single- 
point upwind difference scheme modelling the derivative; i.e., 



/d<M SPUDS 2 K i + H ” 6 < t > i - 1 + *1-2 

te), 6 h 

(42) 

This is then written in /wewdo-flux-difference form 



/ d<t> \ spuos _ (*;) spuds - (*;) spuds 
\di)i h 

(43) 

where 

«.; ) spra ’ s = !«■,., ♦ *,) - ±»,., - 2*1 - 

(44) 

and 

(*;) s,uds = i(*. * - !«■. - 2*,., * 

(45) 


Johnson and MacKinnon call this a “finite-volume” formulation 4 . When this convection 
model is used in Equation (37), using the exact diffusive flux of Equation (34), the node- 
point error asymptotes to a second - order trend, as seen in the SPUDS(l/6) column of Table 
I. Clearly, SPUDS represents a second - order finit e-volume formulation of the convective 
term. 


Table I. Grid-refinement study of the model convection-diffusion 
equation, using a finite-volume flux-difference OA formulation, 
with exact diffusive fluxes calculated from the analytical solution. 
Two convection schemes are compared. Values shown are those 
of the node-point error at x = 0.75. 


h~ l 

QUICK(l/8) 

SPUDS(l/6) 

4 

-3.93435739 x 10' 2 

-1.33129698 x 10' 2 

8 

-9.24652759 x lO^ 4 

+3.72661451 x lO" 4 

16 

-2.75420514 X 10 -4 

+4.77048555 x 10" 4 

32 

-4.12098809 x 10’ 5 

+ 1.76208352 x 10- 4 

64 

-5.30547084 x 10‘ 6 

+5.08534128 x 10' 5 

RATE 

0(h 3 ) 

0(h 2 ) 
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However, SPUDS was not really designed to simulate ( <t> r - 4>t)/h. Instead, it is a 
model of the single-point differential formulation given by Equation (20), which is then 
written in pseudo-flux-difference form. Testing the scheme with the exact finite-volume 
form of the diffusive fluxes is, therefore, inappropriate. A proper test would use the exact 
second derivative of Equation (26) in pseudo-flux-differenct form. This can be achieved by 
writing 

(d^\ _ (»,')• - (»,T (46) 

\ dx 2 /« h 


where (for the particular model problem under consideration) the right pseudo-gradient is 


(*,V 


(P6 e”* e P ‘ /2 \ ( *V 2 \ 
\ e p£ - 1 / lsinhP A /2/ 


(47) 


and the left pseudo-gradient is obtained by replacing x, by (x, - h). When these formulae are 
used in the hybrid flux formulation, the node-point error shows a third-order trend, as seen 
in the SPUDS(l/6) column of Table II. 

As a matter of interest, using the QUICK(l/8) convection scheme in combination with 
the single-point exact diffusion operator gives an 0(h 2 ) trend. This is shown in the 
QUICK(l/8) column of Table II. This is to be expected from Equation (12), since the finite- 
volume QUICK(l/8) scheme is now being used out of context in a single-point formulation, 
just as the SPUDS(l/6) scheme is 0(h 2 ) when used in a finite-volume formulation. 

Of the two third-order convection schemes — the QUICK(l/8) finite-volume 
formulation in Table I and the SPUDS (1/6) formulation in Table II — note that, the finite- 
volume formulation is asymptotically 33 % more accurate. The reason for this is explained 
by a formal discretisation error analysis in the Appendix. 


Table II. Grid-refinement study of the convection-diffusion 
equation, using the SP differential equation expressed in pseudo- 
flux-difference form, with exact diffusion terms calculated from the 
analytical solution. Two convection schemes are compared. 
Values shown are those of the node-point error at x = 0.75. 


h' 1 

SPUDS (1/6) 

QUICK(l/8) 

4 

-4.35438228 x I0' 2 

-7.45209055 x 10' 2 

8 

-1.86968335 x 10‘ 3 

-2.38132421 x 10' 3 

16 

-3.97474350 x 10- 4 

-1.07787069 x 10' 3 

32 

-5.50611000 x 10’ 5 

-2.71242977 X lO* 4 

64 

-7.06282736 x 10’ 6 

-6.32202472 x 10' 5 

RATE 

OQ?) 

0(h 2 ) 
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Diffusion models 

For evaluating numerical models of diffusion terms, the hybrid finite-volume flux is 

written 


F™(x.) 


I 

h 


e*' n 

Tt 


(e™ - 1) 


- 1) ) model 

" Pi . 


(48) 


with a corresponding formula for = F™ B (x.-h ) . If a single-point formulation 

is used, written in pjendo-flux-difference form, then the appropriate formula (for the current 
exponential solution) is 


F r *(x,) 


h 




P ix, _P./2 


- 1) 


i+X 


model 


(e w - 1) 2 sinh P A /2 


P6 


(49) 


QUICK diffusion flux 

Interpolating a parabola through node values </*,_„ and <£ i+1 , on a uniform grid, 
leads to 

^ /^QUICK _ ~ 4*1 (50) 

Because of a geometric property of the parabola, this is indistinguishable from linear 
interpolation between node values <£, and When used with a finite-volume formulation 
of the exact convection terms, Equation (48), this leads to the 0(rf) convergence shown in 
the QUICK column of Table III. 


Classical central differencing 

If the second derivative at point i is approximated by the second central difference 


d 2 <j> \ modcl _ 4> itl - 2 4 > t + 
dx 2 h h 2 


(51) 


the corresponding pseudo-flux-difference formulation involves the pseudo-g radient 

^ . _ *t*i ~ 

Wr ) - U 


(52) 


which, of course, is superficially identical to the QUICK formulation of Equation (50). But, 
in this case, the appropriate (pseudo) convective flux is given by Equation (49). The grid 
convergence behaviour is, therefore, not identical to that of the QUICK diffusive formula- 
tion. As seen in the CDS column of Table III, single-point classical second-order central 
differencing for diffusion asymptotically generates errors exactly twice as large as the 
identical operator used in a finite-volume formulation, using exact convective terms in each 
case. The reason for this is seen in the Appendix. 
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Table HI. Grid -refinement study of diffusion schemes, using 
exact convection terms. Node-point errors are shown at x = 0.75. 
In Column I, the QUICK finite-volume OA scheme is used. 
Column II gives results for classical central differencing using an 
SP formulation. 


hr 1 

QUICK 

CDS 

4 

+ 1.59466830 x 10' 2 

+3.12478024 X 10' 2 

8 

+4.07354976 x 10' 3 

+8.10497409 x 10‘ 3 

16 

+ 1.02395148 x 10' 3 

+2.04524126 x 10‘ 3 

32 

+2.56328047 x 10- 4 

+5.12489289 X IQ - 4 

64 

+ 6.38006991 x 10' 5 

+ 1.27591015 X 10' 4 

RATE 

0(h 2 ) 

OQi 2 ) 


Fourth-order diffusion terms 

Consider the diffusion operators given by Equations (16) and (17). Table IV shows 
the finite-volume formulation using exact convective fluxes from Equation (48). As 
expected, Equation (17) shows an OQf) trend, whereas Equation (16) — being used out of 
context — is only 0(h 2 ), according to Equation (12). 

Alternatively, Table V shows the single-point formulation using the same diffusion 
operators together with Equation (49) for convection. In this case, the convergence is 
reversed, as expected, since now Equation (17) is being used out of context. 

Note that the fourth-order single-point diffusion operator in Table V generates errors 
more than twice as large as those of the fourth-order finite-volume operator in Table IV. 
This, again, appears to suggest that a true finite-volume formulation is likely to be more 
accurate than the corresponding single-point scheme of the same formal order of accuracy . 


Table TV. Grid -refinement study of diffusion schemes, using a 
finite-volume flux-difference OA formulation, with exact convective 
fluxes calculated from the analytical solution. Two diffusion 
schemes are compared. Values shown are those of the node-point 
error at x = 0.75. 


h' 1 

Equation (17) 

Equation (16) 

4 

-1.74276854 X 10’ 3 

-1.70401744 X 10' 2 

8 

-1.12027932 x 10' 4 

-4.02705160 x 10' 3 

16 

-7.10364599 X 10' 6 

-1.00591394 x 10' 3 

32 

-4.47274921 X 10' 7 

-2.53235117 x 10‘ 4 

64 

-2.79314528 x 10' 8 

-6.33685964 x 10‘ 5 

RATE 

0(h*) 

OQi 2 ) 
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Table V. Grid-refinement study of diffusion schemes, using the 
SP differential equation expressed in pseudo- flux-difference form, 
with exact convection terms calculated from the analytical solution. 
Two diffusion schemes are compared. Values shown are those of 
the node-point error at x = 0.75. 


h' 1 

Equation (16) 

Equation (17) 

4 

-3.81023176 x 10' 3 

+ 1.24370324 X 10’ 2 

8 

-2.56261277 x W* 

+3.78389793 x 10’ 3 

16 

-1.65684268 x 10' 5 

+9.97608855 x 10- 4 

32 

-1.05202438 x 10' 6 

+2.53650099 x 10* 4 

64 

-6.59600528 x 10'® 

+6.35095680 x 10* 5 

RATE 

Olff) 

0(h 2 ) 


Convection-diffusion schemes 

For reference, Table VI shows a grid-refinement study of two convection-diffusion 
schemes: 

(i) The standard finite-volume QUICK scheme, using Equation (6) for convection 
and Equation (50) for diffusion; and 

(ii) The SPUDS scheme, Equation (44), for convection, together with classical 
second-order central-differencing, Equation (52), for diffusion. 

Both schemes are 0(h 2 ) because of the dominance of the diffusion terms at the fine-grid end 
of the spectrum. Note, however, that the SPUDS +CDS scheme asymptotically generates 
errors twice as large as those of the standard finite- volume QUICK formulation. This is 
clarified in the Appendix. 


Table VI. Grid-refinement study of the convection-diffusion 
equation, using a finite- volume flux-difference (or pseudo- flux- 
difference) formulation of two schemes: (i) The standard QUICK 
convection-diffusion scheme, (ii) The SPUDS convection operator 
together with CDS for diffusion. Values shown are those of the 
node-point error at* = 0.75. 


ft' 1 

STANDARD QUICK 

SPUDS + CDS 

4 

+9.96622011 x W 4 

+7.58488389 x 10' 3 

8 

+ 1.25198293 x 10' 3 

+3.87320949 x 10' 3 

16 

+5.56979450 X 10- 4 

+ 1.30280308 x 10' 3 

32 

+ 1.74361676 x 10* 4 

+3.73713298 x 10” 4 

64 

+4.82503616 x 10* 5 

+9.97187451 x 10' 5 

RATE 

<9(ft 2 ) 

0(h 2 ) 












13 


CELL AVERAGES AS DEPENDENT VARIABLES 

Up until this point in the discussion, modelled derivatives and fluxes have been written in 
terms of nodal point values of the dependent variable. But many CFD schemes — especially 
finite-volume formulations — treat the cell averages as .the dependent variables, writing the 
convective and diffusive fluxes directly in terms of the 0’s. [It is also theoretically possible 
to write single-point formulations in terms of cell-averages, although this does not appear 
to have ever been proposed in the CFD literature.] Note that_the distinction only occurs at 
third order and above; for first- and second-order schemes, 4> i = 4> r 

From the analytical solution to the model problem being studied, the exact formula 
for the cell average, defined in Equation (13), can be found as 

- _ e Fix ' [(sinh P A /2)/(P A /2)] - 1 (53) 

e M - 1 

The only difference between this and the exact nodal point solution for Equation (22), 
is the appearance of the hyperbolic-sine factor (in square_brackets). For a given P6, this 
factor depends on the grid size; Table VII shows </>,■ and <f> ( values for h = 1/8. Note that 
the sinh-factor has a Taylor expansion given by 

sinh P a /2 _ , P , l ^ P 1 ^ (54) 

P /2 24 1920 

A 

which should be compared with Equation (13), for example. 


Table VII. Comparison of exact node-point values, </>„ with exact 
cell average values, <t > ( , for the model convection-diffusion problem 
with h = 1/8. 



4>i 


0 

0.125 

0.25 

0.375 

0.5 

0.625 

0.75 

0.875 

1 

0.0 

0.012103427 

0.032058603 

0.064959128 

0.119202922 

0.208635820 

0.356085740 

0.599189560 

1.0 

0.000194956 

0.012424854 

0.032588548 

0.065832859 

0.120643461 

0.211010867 

0.360001531 

0.605645608 

1.010644223 
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Sub-cell interpolation _ 

Given a set of <£ (cell-average) values, one would like to interpolate a sub-cell 4>(x) 
satisfying Equation (13). This can be done quite easily by introducing a discrete integral 
variable, defined by 

+ * i, ( 55 > 


where it is understood that values occur at the right face of the corresponding cell i. An 
interpolation, ip(x), collocated at & values can then be differentiated to give <f>(x): 


4>(x) 


dm 

dx 


(56) 


automatically satisfying Equation (13) — no matter what type of interpolation is used for 
\{/(x), provided the collocation conditions 

\l/(x.+h!2) = for all i (57) 

are satisfied, where is located at the centre of cell i. 

For example, piece-wise polynomial interpolation can be used for interpolating \p{x) 
over each cell i. In this case, discontinuities in slope occur in \j/(x) at cell faces; these 
correspond to discontinuities in value in 4>(x) across cell faces. 

Third-order finite-volume convection terms 

Just as with the QUICK(l/8) finite-volume convection scheme using nodal point 
values, the corresponding third-order cell-average formula is based on local (piece-wise) 
quadratic interpolation of 4>(x). This requires piece- wise cubic interpolation of \p(x). For 
cell i, a symmetrically located cubic is interpolated by collocation through ^_ 2 , ^ f _j, and 
&+i, giving 

m = *, * - * M ) (|) 

- r-- 3 n 3 ^ )( iF 


valid over the range: -h < £ < 0, where £ = x - ( x,+ft/2 ). From its definition, Equation 
(56), the sub-cell interpolation across cell i is then given by the quadratic 

m = \(i M * i t ) ~ |(5 M - 2 i t * i M ) 

+ (i M - it) (|) + jdi. 1 - 2 *i + K l) (|) 


(59) 
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For steady flow (with P6 > 0), the right face value of cell i is given by 

<t> r (i) = <K 1=0) = I(i M + i.) - + ii-i) 


(60) 


This “1/6” formula in terms of cell averages is consistent with the “1/8” formula using 
nodal point values. This can be seen by using Equation (59) to evaluate <f> { in terms of the 


</>’ s: 


*, = <K( = -hl2> = *, - i(* M - 2*, * * M ) 


(61) 


and noting that 1/6 = 1/8 + 1/24. 


Grid-refinement using cell-averages 

If Equation (60) for 4> r (i) and the corresponding formula for 0, given by the 
conservation condition, 4>t(i) = 4> r (z-l), are used together with exact diffusive fluxes, using 
Equations (34) and (15), solution of the resulting difference equation will, of course, 
generate approximate <f> { values (as opposed to <t> ( values). These should be compared with 
the corresponding exact </>,. values by introducing the cell-average error 

CAE,. = 4> i (computed) - (exact) (62) 

using Equation (53) for <j> { (exact) . Table VIII shows CAE values atx = 0.75 for convective 
modelling errors using Equation (60), together with exact diffusive fluxes. As expected from 
the quadratic construction of face values, this is a third-order accurate scheme. One can, 
of course, retrieve node values by using Equation (61). The corresponding NPE atx = 0.75 
is also shown in the table. 


Table Vm. Grid refinement study of the finite-volume convection 
scheme corresponding to Equation (60), based on cell-average 
values and using exact diffusive fluxes . The cell-average error is 
shown together with the corresponding node-point error using 
Equation (61). 


h- 1 

CAE (@ x = 0.75) 

NPE (@ x = 0.75) 

4 

-4.53809633 x 10* 2 

-5.18134013 x 10' 2 

8 

-1.88922017 x 10' 3 

-1.90618632 x 10' 2 

16 

-3.98510249 x 10- 4 

-4.48060981 x 10' 4 

32 

-5.50969540 x 10' 5 

-6.76253568 x 10‘ 5 

64 

-7.06397697 x 10‘ 6 

-1.01264031 x 10' 5 

RATE 

0(h 3 ) 

0(6’) 
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Single-point formulation using cell-averages 

For completeness, the single-point formulation using cell-averages should be con- 
sidered. To achieve third-order accuracy, this requires an upwind-weighted cubic sub-cell 
reconstruction of <t>(x ) corresponding to a quartic piece-wise polynomial interpolation of 
\p(x), collocated at ^,_ 3 , &_ 2 , and (for Pd > 0). This gives 

/tf<M modcl _ + 15 ?« " 27 ^»-i + 5 ^»- 2 

12F/, 24/1 

This can be expressed in pseudo- flux-difference form by identifying 

= |(5w + 5|) " ^ (J M ~ 2 ^< + 

with 4>t(i) = 4>' O'-l) , as usual. Table IX shows the corresponding grid convergence of 
CAE and NPE at x = 0.75. 


(63) 

(64) 


Table IX. Grid refinement study of the SP convection scheme 
corresponding to Equation (64), based on cell average values and 
using exact diffusion terms. The cell-average error is shown 
together with the corresponding node-point error using Equation 
(61). 


/ 

CAE (@x = 0.75) 

NPE (@x = 0.75) 

4 

-4.86965236 x 10’ 2 

-5.52922960 x 10' 2 

8 

-2.86590535 x 10‘ 3 

-2.90394004 x 10‘ 3 

16 

-5.11359278 x 10- 4 

-5.62692325 x 10- 4 ! 

32 

-6.86685845 x 10' 5 

-8.11969093 x 10' 5 

64 

-8.81228514 x W 6 

-1.18744119 x 10- 5 

RATE 

0{h 3 ) 

0(h 3 ) 
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CONCLUSION 

In constructing convection-diffusion schemes, there are four general categories. First, one 
can choose to model first and second derivatives at a single point; this is the SP formulation. 
Alternatively, one can choose to model the face values and gradients in a finite-volume 
formulation; this is the operator-average, or OA, formulation. Each of these formulations 
can use either nodal point values or cell-average values of the dependent variable. Control- 
volume formulations are automatically in conservative flux-difference form. The two SP 
formulations can also usually be written in a psewJo-flux-difference form; in this case, 
conservation is satisfied even through the modelled “fluxes” do not represent the true 
physical fluxes. When working with third (or higher) order convection methods, it is 
important to model all terms in the equation in the same manner (i.e. , either all SP or all 
OA); a mixed formulation is condemned to (at best) second order accuracy, no matter how 
accurate the individual terms. 

For a model problem with a known analytical solution, the order of accuracy of the 
convection terms in isolation can be studied in a grid refinement test, using exact diffusion 
terms (and vice versa). Tests of this type showed that, for methods using nodal point values: 

• QUICK(l/8) is 0(h 3 ) accurate in an OA formulation but only 0(h 2 ) 
accurate in an SP formulation. 

• SPUDS(l/6) is 0(rf) accurate in an SP formulation but only 0(h 2 ) 
accurate in an OA formulation. 

For methods using cell-average values as dependent variables: 

• The “1/6” formula gives 0(?P) convergence in an OA formulation. 

• An SP formulation requires a “5/24” factor for 0(rf) accuracy. 

In general, finite-volume formulations are considerably more accurate than the 
corresponding SP formulation of the same formal order. In particular, the QUICK(l/8) 
convection-diffusion scheme in a finite-volume formulation is asymptotically twice as 
accurate as using SPUDS(l/6) for convection with CDS for diffusion in an single-point 
formulation, although both schemes are formally only second-order accurate because of the 
dominance of the diffusion terms in the fine-grid limit. 
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APPENDIX 

RELATIONSHIP BETWEEN TRUNCATION ERROR AND 
DISCRETISATION ERROR 

In the following analysis, it is assumed that the problem is linear, steady, and one- 
dimensional, involving a uniform spatial grid of mesh-size h. Generalisation of these 
conditions are relatively straight-forward. The analysis establishes the (often taken for 
granted) fact that discretisation error in the solution is of the same order of accuracy as 
truncation error in the operator, in a single-point formulation. The same is true of finite- 
volume formulations. 

Single-Point Formulation 

Consider a linear operator involving derivatives of various orders at a single point. This 
(exact) differential operator is represented by L. Assume that U(x) is the exact , 0(1), 
solution of the exact differential equation 

L(U) = S(x ) (A,1) 

where S(x) is a known source term. Now consider a discrete operator representing a 
numerical approximation to the single-point differential operator; this is represented by D. 
Assume that «, is the exact nodal-point solution of the approximate numerical difference 
equation, defined at nodal points i; i.e., 

D(u t ) = S { ( A - 2 > 

where S t = S(x t ) is known exactly (an approximate 5, will, of course, introduce “source- 
term” errors). 

To define the truncation error of the numerical operator, assume that v(x) is a test 
function with all derivatives. Then the truncation error can be defined as the difference 
between the approximate and exact operators, operating on v at the nodal point i: 


TE sp (v) | f = D(v) | . - L(v) | . 


(A. 3) 


By making Taylor series expansions of the terms in D(v)| f , the truncation error will be 
found to depend on h, as follows 

TE sp (v ) | . = P(v)|,. h p + HOT = 0(h p ) (A- 4 ) 

where P(v) involves derivatives of v, and p is an integer. Note that using U as the test 
function in Equation (A. 3) gives 

D(V) |, = 1(17)|, - TE sp (t7) |, = S, *TE sp (t7)|, (A.5) 

using Equation (A. 1), evaluated at i. This equation is valid only at points where U and all 
its derivatives are continuous. Singular points require special treatment, and are excluded 
from the following analysis. 
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Define the nodal-point discretisation error as the difference between the exact nodal- 
point solution of the approximate numerical equation and the exact solution of the exact 
differential equation, evaluated at i: 



(A. 6) 


This is the same as the node-point error, defined previously. [Note that the present analysis 
uses a consistent definition of ‘ ‘error’ ’ as 

“error” = “approximate” - “exact” (A. 7) 

Some authors sometimes use the reverse (negative) of this, which strictly should be called 
“correction’ ’ rather than ‘ ‘error’ ’ .] 

Consider the discrete operator applied to the discretisation error (using the assumed 
linearity property) 

D(e,) = D(u,) - D(U,) < A - 8 > 


or, from Equations (A.2)-(A.5), e t satisfies the discrete equation 

Die,) = -TEgpflOl, = -P(U)\, h p + HOT (A-9) 

Note that S, from Equations (A.2) and (A.5) has cancelled; otherwise, additional “source- 
term” errors must be retained. Now rewrite the left-hand side as 

Die,) = £(«)!, * [B(«.) - 1(e) |,] = m\, * TE sp (e)|, (A. 10) 


where e(x) is a continuous function, with e. = e(x ,) . This means that e(x) satisfies the 
differential equation 


L{e)\ K = - TE SP (l/) 1 1 - TE SP (e)|. 


(A. 11) 


where i x is the nodal value of a fixed point, x = const, as h is varied, 
leading order, 

e = 0(h«) 


Assume that, to 
(A. 12) 


where q > 1. Then Equation (A. 11) becomes 

L{ e) | . = - p(U ) \.h p + HOT + 0{h p * 9 ) 

= 0(h p ) 


(A. 13) 


And, since L is a linear homogeneous operator, independent of h, this means that the leading 
single-point discretisation error is 


e fJ = 0(h p ) 


(A. 14) 
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i.e., that q = p. The discretisation error of a discrete operator in a single-point formulation 
is thus of the same order as the single-point truncation error. 


Finite-Volume Formulation 

Assume that the exact differential operator given by Equation (A. 1) is averaged over a finite- 
volume cell, i. The corresponding (exact) finite-volume equation is then 

L(U) = S(x ) (A ’ 15) 


where, in general, the operator average is 


L 


1 

VOL 


III 


L d V 


(A. 16) 


or, in one dimension, 


L 


x*hJ 2 

— [ L dx 

h J 

x-k/2 


(A. 17) 


Again, U(x ) is considered to be the exact solution of the exact finite-volume equation, where 
S(x) is the known cell-average source term. 

Let the corresponding approximate numerical finite-volume difference equation be 
represented by 

D{u t ) = S f < A - 18 ) 

where u t is the exact nodal-value solution of this approximate equation, and S t is the known 
exact cell-average source term at cell i. Note that finite-volume (or, for that matter, single- 
point) formulations can be written in terms of cell-average values, w, , rather than nodal 
values, In general, this will involve a different discrete operator. The present analysis 
will focus on nodal values, but entirely similar conclusions can be reached using cell-average 
values. 

Once again, the truncation error of the discrete finite-volume operator is defined as 
the difference between approximate and exact operators, operating on a test-function, v, at 
cell i: 

TE FV (v) | ~ = D(v)\ i -Z(v)| i | (A- 19) 


Taylor series analysis leads to 

TE fv (v)|, = *(v) |,ft' * HOT 

analogous to Equation (4). Using U instead of v gives 

D(U)\. = L(U)\ i + TE fv (I/)|j = S t + TE fv ( 77) | f 


(A. 20) 
(A. 21) 


for each cell, i. 

The nodal-point discretisation error is again defined by Equation (A. 6). This now 
satisfies the following discrete equation 
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D(e,) = -TE FV (l/)|, = ”11(12)1; h r * HOT (A.22) 


Using Equation (A. 10) results in a differential equation for e : 

L(e)\. = -TE FV (U)|. -TE SP (e)|. 


(A. 23) 


Note that the left-hand side has been written in terms of L rather than L , since L depends 
on h. Assume that, to leading order, 


e = 0(h r ) 


(A. 24) 


Then Equation (A. 23) becomes 

L(e ) |,. = -/?(U)| f h r + HOT + OQi r * p ) 

= 0(h r ) 


(A. 25) 


as assumed, since p > 2, according to Equation (12). The discretisation error of a discrete 
operator in a finite-volume formulation is thus of the same order as the finite-volume 
truncation error: 



= 0(h r ) 


(A. 26) 


Note, however, that for a given discrete operator, treated alternatively as a single- 
point or a finite-volume operator, the corresponding respective truncation (and, hence, 
discretisation) errors will be different. In fact, according to Equation (12), they will differ 
by 0(h 2 ). 

Finally, it should be clear that, for any two different discrete operators or different 
(SP or FV) formulations representing a physical quantity (e.g., convection), the ratio of the 
discretisation errors will be in the same proportion as the ratio of the respective truncation 
errors, as h -* 0. This will be demonstrated in the following section. 

EXAMPLES 

The model convection-diffusion problem introduced earlier forms a good example for 
studying the relationship between truncation error and discretisation error. For convenience, 
the governing equations are repeated here; written in terms of the exact solution, U(x). 

SINGLE-POINT FORMULATION: 


dU _ 1 d 2 U 
dx P6 dx 2 


0 


(A. 27) 
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FINITE- VOLUME FORMULATION: 



(A. 28) 


BOUNDARY CONDITIONS: 


U( 0 ) - 0 , 


17(1) = 1 


(A. 29) 


EXACT SOLUTION: 


U(x) = 


.Viz 


- 1 


e M - 1 


(A. 30) 


DERIVATIVES: 


dU 

dx 


P 6 e 




- 1 


d 2 U 

dx 2 


V€ e 


Viz 


,V 6 - 1 


(A. 31) 


FUNDAMENTAL FV-SP RELATIONSHIP: 



(A. 32) 


Convection Operators 

The single-point upwind difference scheme for convection recommended in Reference 4 can 
be written 


SPUDS 


2 4>.- +1 + 3< £, " 6< £i-l + 4>t- 2 

6h 


(A. 33) 


In terms of a test-function, v, this has a Taylor expansion about grid-point i as follows 


SPUDS 


dv 

dxh 


V . (iv) v (v) v (v!) 

+ _L_ h 3 - — h* + 

12 30 72 


V . (vii) 

h s - — h 6 

252 


(A.34) 


This is the appropriate form for a single-point formulation. Using Equation (A. 32), for 
n = 0, gives the corresponding finite- volume formulation 


SPUDS 





I3v/ V) 

nr 


7t 4 


+ 



h 5 


61v/ vii) 

15360 


(A. 35) 



23 


The QUICK convection operator can be obtained from Equations (40) and (41) as 


■ QUICK 
<Pr 


■ QUICK 
“ <Pt 


[qUICK(C)] = 

The Taylor expansion about grid-point i is 

[ QUICK(C) 


3<*> i+1 + 3<f>, - 7^ f 4>,. 2 (A. 36) 


8/t 


*Z) + Zi_ 

dxli 24 


v (iv) llv' 

*’*Tr*’- 


(V) 


v (v0 

+ — — h* 

480 96 


59v/ v “ > 

20l60 


(A. 37) 


This would be the form used in a single-point formulation. However, QUICK is specifically 
designed for a finite-volume formulation; using Equation (A.32), for n = 0, leads to 


[ QUICK(C) ] 



V: 


(iv) 


3v (v) v- 

+ — h s 
128 96 


(Vi) 


16 


h 3 - 


3v/ vii) 

1024 


h< 


(A. 38) 


Numerical Values 

To get some idea of the relative size of truncation error terms, the known exact solution of 
the model problem is used, with h = 1/64. This would normally be considered a “very 
fine” grid; and asymptotic trends have been established, as seen in the previous tables. The 
truncation error for the single-point formulation of the SPUDS operator is derived from 
Equation (A. 34) using U in place of v: 


jj 9 v) 

TE sp [SPUDS] = -L _ - _1_ i 4 * A 5 


US 


(vii) 


252 


h ( 


(A. 39) 


For the particular model problem being considered, this gives, at x = 0.75, 


TE sp [SPUDS] 


(2£ h> - h* 
\ 12 30 



P6 7 

252 


h 



e 015 ?i \ 
e™ - II 


(A. 40) 


For Pe = 4 and h — 1/64, the numerical values of the individual terms are, respectively 
TE sp [SPUDS] = (8.13802 - 0.20345 + 0.00530 - 0.00010) 3.74743 x 10’ 6 (A.41) 


or 

TE sp [SPUDS] = 2.97538 x 10' 5 (A.42) 

Note that the second, 0(h A ), term in the truncation error is not insignificant. 

If SPUDS is now considered (albeit inappropriately) as a finite-volume formulation, 
Equation (A. 35) gives the corresponding truncation error as 
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TE fv [SPUDS] = - 


JJ"t TjO") 

Ui h 2 t -L- h 3 - 


24 ' 12 

ewr 


13U? fc4 ^ V f* 

384 + 72 


15360 


h 6 + ... 


(A. 43) 


with numerical values 

TE fv [SPUDS] = (- 65.10417 + 8.13802 - 0.20663 + 0.00530 - 0.00002) 

x 3.74743 X 10' 6 


(A. 44) 


or 


TE pv [SPUDS] = - 2.14231X10- 


(A.45) 


From Equation (A. 37), the QUICK convection scheme used (inappropriately) in a 
single-point formulation would have a truncation error 


V" „ U[ n) , 11 Ui 

TE SP [QUICK(C)] = -J— h 2 + — ^ k 3 - 


24 


16 


480 


(V) tt(v0 

h 4 + 


96 


r<vii) 


59U> 

l6 


20160 

giving numerical values 

TE SP [QUICK(C)] = (65.10417 + 6.10352 - 0.13987 + 0.00397 - 0.00007) 

x 3.74743 x 10* 6 


(A. 46) 


(A. 47) 


or 


TE SP [QUICK(C)] = 2.66336 x HT 4 
By contrast, the (appropriate) finite-volume formulation of QUICK leads to 


TE PV [QUICK(C)] = 


u[ iv) 3 1 /, (v) , 4 t// vi) , , 3l// vii) 


16 


hr - 


128 


ft 4 + 


96 


ft 5 - 


1024 




(A. 48) 


(A. 49) 


with numerical values 

TE fv [QUICK(C)] = (6.10352 - 0.14305 + 0.00397 - 0.00007) 

x 3.74745 X 10- 6 

or 

TE PV [QUICK(C)] = 2.23511 X 10‘ 5 


(A. 50) 
(A. 51) 
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Recall from Tables I and II, node-point errors for h — 1/64: 

NPE fv [QUICK(C)] = - 5.30547084 x 10' 6 ( A - 52 ) 

NPE pv [SPUDS] = + 5.08534128 x 10' 5 (A.53) 

NPE sp [SPUDS] = - 7.06282736 X 10’ 6 (A.54) 

NPE sp [QUICK(C)] = - 6.32202472 x 10‘ 5 ( A - 55 ) 

Note, in particular, the ratio of the single-point SPUDS error to the finite-volume QUICK(C) 
error (i.e., the two third-order methods): 

NPE sp [SPUDS] _ 133 (A.56) 

NPE fv [QUICK(C)] 


This is virtually the same as the ratio of the respective leading truncation error terms. From 
Equations (A. 39) and (A. 49), 

LTE sp [SPUDS] = 4 (A51) 

LTE fv [QUICK(C)] 3 

This relationship will be found to be (approximately) true for other ratios, as well. 


Diffusion Operators 

The second-order central-difference operator for diffusion is 


CDS2 


1 _ 2 4>, + ^«-i 

h 2 


(A. 58) 


This can be viewed as a single-point operator 


CDS2 


] - m 


v; 


(iv) 


v; 


(v0 


V (viii) 

1,6 


12 360 20160 
or a finite-volume operator, using Equation (12), with n — 1, 


CDS2 


v/ - v/ 


v f (W) . „ I3v/ Vi) 


24 


h 2 + 


Vi 


(viii) 


5760 


21504 


(A. 59) 


(A. 60) 


Note that the latter form (with the smaller truncation error) corresponds to the QUICK finite- 
volume formulation of diffusion, QUICK(D). 

For Pe = 4 and h = 1/64, the respective numerical values are 

TE sp [CDS2] = (5.20833 x 10‘ 3 + 6.78168 X 10’ 7 + 4.73053 x 10' 11 ) 

x 3.74743 x IQ' 1 (A ‘ 61 ^ 


and 
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TE fv [QUICK(D)] = (2.60417 X 10’ 3 + 5.51012 x 10' 7 + 4.43487 x 10' 11 ) 

x 3.74743 xlO' 1 

In either case, the leading term is dominant. Note that 

LTE sp [CDS2] _ 0 
LTE fv [QUICK(D)] 


This is reflected in the node-point error of Table III, where 

NPE fv [QUICK(D)] = + 6.38006991 x 10’ 5 


(A. 62) 
(A. 63) 
(A. 64) 


and 

NPE sp [CDS2] = + 1.27591015 x KT 4 ( A - 65 ) 


This means that the finite-volume (QUICK) formulation of diffusion is twice as accurate as 
the single-point CDS2 formulation, although both are formally second-order accurate. 

The fourth-order single-point diffusion operator given by Equation (16) has a single- 
point truncation error 


TE sp [Eq(16)] 



t7/ vii ° 


1008 


(A. 66) 


When used (inappropriately) as a finite-volume operator, this gives 


TE FV [Eq(16)] 


u; 


(iv) 


24 


h 2 - 


13[/, <vi) 

1152 


h 4 


10717/™° , , 

h 6 + 

107520 


(A. 67) 


which, of course, is only second-order accurate, according to Equation (12). 

By contrast, the fourth-order finite-volume operator given by Equation (17) has a 
finite-volume truncation error 


TE FV [Eq(17)] 


3U/ Vi) 

“640" 


h 4 


19317/ “° ,, 

: — h 6 + 

322560 


(A. 68) 


And if this were used (inappropriately) in a single-point formulation, the truncation error 
would be 


TE SP [Eq(17)] 



u; 


(Vi) 


240 


1917/™° 

1080 


h 6 


+ . . . 


(A. 69) 


again, second-order, according to Equation (12). Note that 

LTE sp [Eq(16)] _ ^ 
LTE FV [Eq(17)] 


(A. 70) 
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From Tables IV and V, for h = 1164, 

NPE sp [Eq(16)] = - 6.59600528 X 1Q-* = 7 16 (A.71) 

NPE pv [Eq(17)] - 2.79314528 x 10' 8 

Again, one sees that the finite- volume fourth-order formulation is significantly more accurate 
than the single-point fourth-order formulation. 


Convection-Diffusion Operators 

The standard QUICK scheme for both convection and diffusion can be written 


[quick] = 1 

[3* w + 3 <f>i - 7* M + 

j- ±1 

* 2< t>i + 4>t- 1\ 

i 8 h 1 

1 Pe 1 

i h 2 ' 


Viewed as a finite-volume formulation, this gives 

QUI ck] = 


1 /v,- 


(iv) 


Pe\ 24 

h 


h 2 + 


v, 


(iv) 


v (vi \ 

v > \ ij 


96 


16 

3v/ vii) 


h 2 - 
1 


128 


(viii) 


v« * J_( 13V ' 


(vi) 


Pe \ 5760 


1 1024 " Pe I 21504 ). 


h 6 + 


Whereas a single-point analysis would give 

[<=«*]- a - r(£), 


v"' 1 v (iv) \ 

- J_ _L_U 2 


24 


Pe 12 


(iv) 

V. \ 

ir 1 *’- 


iiv; 


(V) 


1 


L 480 Pe V 360 


,(* 0 ,, 


v, 


(vi) 


96 


h 5 - 


59v/ vi0 + j_ I vf* V 


20160 P6 \20160/. 


(A. 72) 


(A. 73) 


(A. 74) 


Although the truncation error terms look different, it is not hard to show that, in fact, in 
terms of U, they are identical. For example, since U satisfies the exact equation 


then 


dU = J_ d^U (A. 75) 

dx Pe dx 2 

jj>" — 1 £/(* v ) (A. 76) 

Pe 
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and the leading single-point truncation error term becomes 


M'" _ i 

U (iv) \ 

• U 2 = 

(2. 

u/ iv) 
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l 24 Pe 

12 / 

\P6 

24 

Pe 

12 / 



and, similarly for the other terms. 

The SPUDS-plus-CDS2 convection-diffusion scheme can be written 

[ SPUDS +CDS2 ] - ( 2 *,. - £ ( *'" - 

As a single-point formulation, this gives 


1 — l^ v \ 

1 

tdh\ 

J \dxli 

Pe 

i& 2 <. 
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whereas, for a finite-volume formulation, 


[ SPUDS +CDS2 | = 


v. - v. 


J_/v/ - v/ 
P6 


24 
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(A. 78) 


(A. 79) 
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(A. 80) 


Once again, in terms of U, the two truncation errors are identical. 

Note that, because of the dominance of the second-order diffusion terms, as h -* 0, 

LTE[SPUpS +CDS2] = 9 (A. 81) 

LTE[QUICK] 

This is borne out in Table VI, where it is seen that, for h = 1/64, 

NPE[SPUDS +CDS2] = 9.971 87451 X 10~ s = 2 07 (A. 82) 

NPE [QUICK] 4.82503616 x 10' 5 
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Finally, it should be noted that in high-convection problems, where the grid Peclet number 

P. = h Pd ( A - 83 ) 

A 

is large, the appropriate way to write the leading truncation error terms is as follows: 

Tj(* v ) i i 

LTE PV [QUICK] = -JL_(l - ^Ljh 3 ( A -84) 

and 

LTE sp [SPUDS +CDS2] = _^_|l - J_p 3 ( A - 85 ) 

This means that, for most flows of interest, using practical grids (so that P A is very large), 
these convection-diffusion schemes are effectively third-order accurate (even though only 
formally second-order accurate as h -* 0). 



30 


REFERENCES 

1. B.P. Leonard, ‘A stable and accurate convective modelling procedure based on 
quadratic upstream interpolation’, Comput. Methods Appl. Mech. Eng., 19, 59-98 
(1979). 

2. B.P. Leonard and S. Mokhtari, ‘Beyond first-order upwinding: the ULTRA-SHARP 
alternative for nonoscillatory steady-state simulation of convection’, Int. J. Numer. 
Methods Eng., 30, 729-766 (1990). 

3. B.P. Leonard, ‘Comparison of truncation error of finite-difference and finite-volume 
formulations of convection terms’, NASA TM 105861, ICOMP-92-19, NASA Lewis 
Research Center (1992). 

4. R.W. Johnson and R.J. MacKinnon, ‘Equivalent versions of the QUICK scheme for 
finite-difference and finite-volume numerical methods’. Comm. Appl. Num. Methods, 
8, 841-847 (1992). 

5. D. Bradley, M. Missaghi and S.B. Chin, ‘A Taylor-series approach to numerical 
accuracy and a third-order scheme for strong convective flows’, Comput. Methods Appl. 
Mech. Eng., 69, 133-151 (1988). 

6. C.A.J. Fletcher, Computational Techniques for Fluid Dynamics, 2nd ed., Vols. I and 
II, Springer, New York (1990). 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden lor this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway. Suite 1204. Arlington. VA 22202*4302. and to the Office of Management and Budget. Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 


November 1993 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


Order of Accuracy of QUICK and Related Convection-Diffusion Schemes 



7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSEES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


5. FUNDING NUMBERS 


WU - 505-90-5K 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-8236 


9. SPONSORING/MO NITORING AGENCY NAME(S) AND ADDRESSEES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSOR1NG/MONITORING 
AGENCY REPORT NUMBER 

NASA TM- 106402 
ICOMP-93-47 


11. SUPPLEMENTARY NOTES 

BP. Leonard, Institute for Computational Mechanics in Propulsion, NASA Lewis Research Center, and The Univer- 
sity of Akron, Center for Computational Mechanics, Akron, Ohio 44325-3903, (work funded under NASA Coopera- 
tive Agreement NCC3-233). ICOMP Program Director, Louis A. Povinelli, (216) 433-5818. 


12a. DISTRIBUHON/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Category 64 


13. ABSTRACT ( Maximum 200 words) 

This report attempts to correct some misunderstandings that have appeared in the literature concerning the order of 
accuracy of the QUICK scheme for steady-state convective modelling. Other related convection-diffusion schemes 
are also considered. The original one-dimensional QUICK scheme written in terms of nodal-point values of the 
convected variable (with a 1/8-factor multiplying the “curvature” term) is indeed a third-order representation of the 
finite-volume formulation of the convection operator average across the control volume, written naturally in flux- 
difference form. An alternative single-point upwind difference scheme (SPUDS) using node values (with a 1/6- factor) 
is a third-order representation of the finite-difference single-point formulation; this can be written in a pseudo-fl ux- 
difference form. These are both third-order convection schemes; however, the QUICK finite- volume convection 
operator is 33% more accurate than the single-point implementation of SPUDS. Another finite-volume scheme, 
writing convective fluxes in terms of cell-average values, requires a 1/6-factor for third-order accuracy. For complete- 
ness, one can also write a single-point formulation of the convective derivative in terms of cell averages, and then 
express this in pseudo-flux-difference form; for third-order accuracy, this requires a curvature factor of 5/24. Diffu- 
sion operators are also considered in both single-point and finite- volume formulations. Finite-volume formulations 
are found to be significantly more accurate. For example, classical second-order central differencing for the second 
derivative is exactly twice as accurate in a finite-volume formulation as it is in single-point. 


14. SUBJECT TERMS 


Quick truncation error; Finite-volume; Operator-average; 
Finite-difference; Single-point 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

32 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Sid. Z39-18 
298-102 

























